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ABSTRACT 


>  this  report,  a  large  number  of  motion  measurements  from  con¬ 
tained  explosions  in  various  rock  types  were  reviewed  and  presented.  These 
data  were  interpreted  in  terms  of  dimensional  analysis  and  other  constraints 
necessary  for  dynamic  consistency  to  yield  scaling  relations  for  estimating 
peak  motions  in  various  rock  media  from  contained  bursts.  Suggestions  are 
also  given  for  estimating  ground  motions  from  direct  induced  ground  shock 
resulting  . rom  a  surface  burst  on  rock. 
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PREFACE 


This  investigation  was  authorized  by  the  Chief  of  Engineers  (DAEN-MCE-D) 
and  was  performed  in  FY  71  and  72  under  contract  No.  DACA  1*5-69-0-0100  between 
the  Omaha  District,  Corps  of  Engineers  and  Dr.  A.  J.  Hendron,  Mahomet,  Illinois. 
The  work  is  part  of  a  continuing  effort  to  develop  methods  which  can  be  used  to 
design  underground  openings  in  jointed  rock  to  survive  the  effects  of  nuclear 
weapons. 

This  report  was  prepared  under  the  supervision  of  Dr.  A.  J.  Hendron,  Princi¬ 
pal  Investigator.  During  the  work  period  covered  by  this  report.  Colonel  B.  P. 
Pendergrass  and  Colonel  Alfred  L.  Griebling  were  District  Engineers;  R.  G.  Burnett 
was  Chief,  Engineering  Division,  C.  J.  Distefano  was  Technical  Monitor  for  the 
Omaha  District  under  the  general  supervision  of  Kendall  C.  F^x,  Chief,  Protective 
Structures  Branch.  Dr.  J.  D.  Smart,  R.  G.  Dodson,  and  D.  G.  Heitmann  particv^ted 
in  the  monitoring  work. 
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CONVERSION  FACTORS ,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 


British  units  of  measurement  used  in  this  report  can  he 
metric  units  is  follows; 


Multiply 

inches 

feet 

cubic  inches 
pounds 

pounds  per  square  inch 


0. 3048 
16.3871 
0.45359237 
0.070307 


pounds  per  cubic  foot 


lt>.  01 85 


inch-pounds  0.0U52I 

inches  per  second  2.54 

inches  pep  second  per  aeeoud  2.54 


vii 


converted  to 

To  Obtain 

centimeters 

meters 

cubic  centimeters 

kilograms 

kilograms  per  square 
centimeter 

kilograms  per  cubic 
meter 

Meter-kilograms 

centimeters  per  second 

centimeters  per  second 
pep  second 


Introduction 


The  design  of  close  in  protective  structures  in  rock  requires  the 
estimation  of  the  ground  motions  from  direct  ground  shock.  At  present,  the 
magnitude  of  these  motions  is  best  estimated  by  scaling  the  field  measure¬ 
ments  from  fully  contained  nuclear  detonations.  The  fully  contained  data 
are  then  adjusted  by  means  of  a  '‘coupling"  factor  to  apply  to  the  surface 
burst  case  considered  for  design. 

In  this  report,  the  concepts  or  scaling  are  discussed  and  field 
roaasuresfients  from  underground  explosions  are  scaled  according  to  the  prin¬ 
ciples  developed.  Recommendations  are  also  given  for  extrapolation  of  the 
contained  data  to  estimate  motions  in  a  rock  medium  below  a  surface  burst. 

Dimensional  Analysis 

The  empirical  scaling  of  shock  phenomena  from  explosions  involves 
the  comparison  of  dynamic  measurements  obtained  at  var^r?  4l&tweft&  in  dif - 

■v 

ferent  media  from  a  wide  range  of  yields.  i]ui te  often  scaling  is  used  to 
estimate  ground  motions  by  extrapolation  to  yields  which  are  far  beyond  the 
range  of  yields  from  which  the  experimental  data  were  obtained.  For  these 
cases,  ft  is  very  important  that  the  scaling  relation  be  guided  by  dimen¬ 
sional  analysis  such  that  the  empirical  relations  obtained  are  dimensionally 
homogeneous.  It  is  also  important  that  Certain  relations  exist  among  empir¬ 
ical  equations  for  maximum  acceleration,  particle  velocity,  and  displacement 
such  that  the  equations  are  kinematically  consistent.  In  this  section,  prin¬ 
ciples  are  developed  for  formulating  empirical  equations  for  ground  motions 
from  experimental  data  which  are  kinematically admissible  and  which  are  di- 
mens  i  cma  lly  homogeneous  . 
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The  dependent  variables  considered  in  such  an  analysis  are  the 
displacement  6,  the  particle  velocity v,  and  the  acceleration  a,  which  result 
at  some  point  due  to  the  explosion.  The  independent  variables  considered 
the  most  significant  in  influencing  the  ground  motions  are  the  energy  re¬ 
leased  by  the  detonation  W,  the  range  from  the  detonation  to  the  point  of 
observation  R,  the  density  of  the  rock  p,  the  compressional  seismic  velocity 
in  the  rock  mass  c,  and  time  t.  Table  I  gives  dimensions  of  all  variables 
considered  in  terms  of  force  F,  length  L,  and  time  T. 

According  to  the  Buckingham  Pi  theorem  (Murphy,  1950),  the  depen¬ 
dent  variables  are  related  to  dimensionless  groups  of  the  independent 
variables  as  follows: 
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whore  g-j ,  g^,  and  g^  are  unknown  functions  of  the  dimensionless  groups  of 
independent  variables  (Pi  terms). 

If  the  maximum  values  of  displacement,  particle  velocity,  and 
acceleration  are  of  interest,  the  Pi  term  tc/R  is  usually  neglected  because 
this  term  relates  only  to  the  time  scaling  of  the  phenomena.  It  is  inter¬ 
esting  to  note  however,  that  R/c  is  the  transit  time  of  a  pulse  from  the  point 
of  detonation  to  a  point  at  a  distance  R.  Thus,  the  Pi  term  tc/R  indicates 
that  "times"  (rise  time,  duration  time,  etc.)  describing  the  character  of 
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-the  pulse  should  scale  directly  proportional  to  the  transit  time.  Figure  1 

shows  a  typical  qualitative  wave  form  produced  by  shock  directly  transmitted 

into  rock  by  explosions.  In  general,  available  data  indicate  that  the  rise 

time  to  peak  particle  velocity  tr,  as  shown  on  Figure  1,  is  1/6  to  1/12  the 

transit  time  of  the  pulse  (Newmark  and  Haltiwanger,  1962).  The  positive 

phase  duration  of  the  velocity  pulse  which  corresponds  to  the  time  of 

maximum  displacement  is  roughly  1  to  2  times  the  transit  time  as  illustrated 

in  Figure  1.  Thus,  the  available  field  data  appear  to  support  that  "times" 

scale  proportional  to  the  transit  times  as  indicated  in  Eqs.  1  through  3. 

Thus,  if  maximum  values  of  6,  v,  and  a  occur  at  constant  values  of  tc/R,  or 

if  tc/R  is  not  important  in  determining  the  maximum  values  of  6,  v,  and  a, 

then  the  dimensionless  maximum  displacement,  velocity,  and  acceleration  are 

?  3 

given  as  functions  g^,  gy,  and  gz  of  only  the  dimensionless  variable  W/pc  R 
as  given  below. 


Thus,  according  to  the  dimensional  analysis  presented,  the  field  measurements 
of  maximum  acceleration,  velocity,  and  displacement  at  various  ranges  from 
different  yield  detonations  in  different  media  should  be  plotted  as  shown  in 
Figure  2.  If  the  scaling  described  works,  then  the  data  should  collapse 
along  a  given  curve  in  the  plots  shown  in  Fig.  2  and  should  enable  the  func¬ 


tions  g^,  gy,  ar.d  ga  to  be  determined. 


Normally,  the  plots  suggested  in  Figure  2  are  made  on  log-log  paper 
and  much  of  the  relationship  of  interest  can  be  approximated  by  a  straight 
line.  Thus,  the  equations  for  <S,  v,  and  a  are  of  the  form: 
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Thus: 


6  =  K  R 


V  lol  V3  c2/3  "6 


V  =  K-|R  v 
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According  to  dimensional  analysis, the  exponents  n*,  n„,  and  n,  in 
Eqs.  7,  8,  and  9  need  not  be  related;  but  since  the  maximum  values  of  dis¬ 
placement,  velocity,  and  acceleration  are  kinematically  related,then  real¬ 
istic  constraints  must  be  investigated  which  are  necessary  to  make  these 
quantities  kinematically  consistent.  For  example,  from  Fig.  1  it  is  apparent 


that  the  maximum  acceleration,  am3„,  is  given  by 

maX 


,  _  »  max 

amax  k3  t 


(13) 


and  that 


l* 


^4  vmax  ud 


where  Kg  and  depend  on  the  shape  of  the  pulse.  For  a  given  shape-  piP-e 


the  rise  time,  t  ,  is  also  given  by 


V  =  S  ^d 


Thus,  substitution  of  Eq.  15  into  Eq.  13  and  multiplication  of  Eqs.  13  and 


14  yields 


a  *6 
max  max 


w2 

Kg  max 


Newmark  (1968)  has  shown  that  the  constant  Kg  K^/Kg  is  always  less  than  0.5. 
Nevertheless,  Eq.  16  shows  that,  independent  of  the  value  of  Kg  K^/Kg, 

Eqs.  7,  8,  and  9  should  satisfy  the  relationship 


2n  =  n,  +  n, 
v  6  a 


in  order  for  Eq.  16  to  be  satisfied. 

For  example,  say  that  it  was  found  that  ny  =  -2  and  na  =  -2  from 
experiments.  For  Eq.  17,  n^  would  also  have  to  be  -2.  Equations  10,  11, 
and  12  then  would  be  of  the  form 

«  .  KR-’f3^3  (18) 

V  =  K.R-2  W2/3  (19) 

'  P 

a  =  K2R-3  f 3  c+’/3  (20) 

4  p 


in  orde>'  to  be  kinematically  consistent. 
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Scaling  of  Field  Measurements  from  Contained  Explosions 


Radial  accelerations  and  particle  velocities  from  the  contained 
explosions  given  in  Table  n  are  given  in  Tables  III,  IV,  and  V.  The  data 
presented  in  Tables  III  -  V  are  shown  scaled  in  the  manner  suggested  in  Fig.  2 
and  plotted  in  Figs.  3,  and  4.  Some  of  the  data  given  are  from  HE  experiments 
and  some  are  from  nuclear  experiments.  It  should  be  noted  that  no  adjust¬ 
ment,  has  been  made  for  the  relative  efficiency  of  nuclear  and  HE  contained 
bursts  for  the  scaled  points  shown  in  Figs.  3  and  4.  If  there  is  a  dif¬ 
ference  in  efficiency  between  HE  and  nuclear  contained  bursts  in  producing 
direct  ground  shock,  then  the  difference  appears  to  be  small  enough  to  be 
masked  by  the  normal  scatter  in  ground  motion  data.  The  relative  position 
of  the  UET  sandstone  data  and  the  Gas  Buggy  sandstone  data  in  Fig.  3  tend 
to  illustrate  the  point  that  there  is  no  discernible  difference  in  efficiency 
which  could  be  detected  within  the  scatter  of  the  data.  From  Figs.  3  and  4, 
the  best  fit  to  the  data  indicates  that  n,  and  n,  =  -2.5.  Thus,  the  form 
of  the  scaling  relation  suggested  is 


then  a  consistent  relationship  for  the  maximum  dynamic  displacement  should 
be  of  the  form: 
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T'ius  a  set  of  consistent  equations  should  be  of  the  form: 
e  -  I  -  W5/6  p'576  c'5/3  *-5/2 

v  =  K^p*5'6  c-2/3  r-5/2 

a  -  K2W5/6p^6c+1/3R'7/2. 

6  -  K  W576  p'576  c"5/3  R“3/2 


(24) 

(25) 

(26) 

(27) 

(28) 


Evaluation  of  the  constants  and  K2  from  the  data  presented  in  Figs.  3 
and  4  yield  equations  for  maximum  radial  strain,  maximum  radial  acceler¬ 
ation  and  maximum  radial  particle  velqcity  as  follows  for-  fully  contained 
nuclear  detonations 
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(32) 


The  maximum  radial  stress,  a  ,  is  given  by 


a  =  pv  c 
r  r  p 


where  cp  is  the  propagation  velocity  of  the  peak  radial  stresses.  In  the 
Climax  Stock  Granite  (location  of  the  Pile  Driver,  Hardhat,  Shoal,  and  Tiny 
Tottests)  the  propagation  of  the  peak  stresses  in  the  close  in  ranges  was 
about  14,000  ft/sec,  and  the  seismic  velocity  was  about  18,000  ft/sec.  Then 
by  £q$.  31  and  32,  the  radial  stress  at  a  range  of  1000  ft  in  a  medium  with 
a  seismic  velocity,  c,  of  18,000  ft/sec  would  be  given  by 


o  =  (2?  ft/sec)(  14,000  ftVsec)  ~  )(  -  )  =  13,500  psi 

r  V:  32.2  ft/sec^  144  lb/fr 

and  a  general  equation  for  radial  stress  would  be  given  by 


a  V;n  -on  nsi  r  w  i5/6r1000  ft-.5/2r  Y _ i1/6r  c _ i1/3 

V  13,l00  psi  ^50  Kt  ^  R  ■*  ^165pcf^  ^18,000  fps^ 


(33) 


where  c  is  the. seismic  velocity  of  the  medium.  Equation  33  is  a  reasonable 
approximation  for  the  radial  stresses  when  the  peak  stress  propagation  velocity, 
cD»  isiabout  (I^qqq  ||^|||  =  0.78)  0.80  times  the  seismic  velocity  of  the 
med  uin.  For  those  cases  where  cp/c  is  different  than  0.80,  the  value  obtained 
from  Eq.  33  should  be  multiplied  by 
c  /c 

or 


C  /c 

(  TpW  )' 


If  the  positive  duration  td  of  the  velocity  pulse  is  approximately 
1  to  2  times  the  transit  time,  R/c,  then  for  1000  ft  from  50  Kt  in  a  medium 
with  a  unit  weight  of  165  pcf  and  a  seismic  velocity  of  18,000  fps  the 
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maximum  displacement  for  a  nearly  triangular  pulse  as  shown  in  Fig.  1  would 
be: 


5max  ~  ^2  vmax 


(34) 


and  K  would  ranye  from 

1/2  x  27  ft/sec  x  18,000  ft/sec  x  1  =  0,75 
to 

1/2  x  27  ft/sec  x  13,000  ft/ sec  x  2  =  1,5  ft 
and  the  equation  for  maximum  radial  displacement  would  be 
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1.5  ft 


j_W_J5/6  1 1000  ftj3/2  1 165  pcf  5/6  |  18,000  ft/secj5/3 


(35) 


for  a  fully  contained  burst. 

Estimation  of  Direct  Induced  Motions  from  Surface  Bursts  on  Rock 


To  estimate  deep  underground  motions,  engineers  have  used  the  con¬ 
cept  of  an  effective  or  equivalent  yield,  W_.  The  equivalent  yield  is  de- 

c 

fined  as  the  yield  of  a  contained  nuclear  explosion  that  would  provide  the 
observed  peak  stress,  strain,  or  particle  velocity  at  a  given  range  beneath 
ground  zero  of  a  nuclear  surface  burst  of  yield  Ws-  The  ratio  Wg/Ws  is  then 
defined  os  an  effective  coupling  factor.  Based  on  peak  particle  velocity 
data  from  low  yield  nuclear  detonations  in  underground  cavities  in  granite 
and  tuff,  an  equivalent  yield  coupling  factor  of  about  10  to  12  percent  would 
be  inferred.  Cooper,  Brode,  and  Leigh  (1971)  have  stated  that  the  "coupling" 
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factors  obtained  from  the  low  yield  cavity  experiments  are  probably  too 
high  to  be  used  for  the  case  of  large  yield  surface  bursts  because  the  more 
massive  (low  yield  to  mass  ratio)  low  yield  devices  are  expected  to  couple 
energy  to  the  ground  more  efficiently  than  the  large  yield  devices.  Thus, 
the  most  current  estimate  of  the  coupling  factor  for  peak  velocities, 
strains,  stresses,  and  accelerations  is  about  5  percent  for  motions  directly 
beneath  a  large  yield  surface  burst  on  a  rock  medium  (Cooper,  Brode,  and 
Leigh,  1971).  An  equivalent  yield  or  "coupling"  factor  however,  would  be 
expected  to  be  lower  for  peak  particle  displacements  than  for  peak  particle 
velocities  because  free  surface  effects  reduce  the  positive  phase  duration 
of  the  particle  velocity  pulse  from  that  which  would  be  observed  in  a  fully- 
tamped  burst.  A  study  of  the  data  for  the  low  yield  cavity  experiments, 
(Cooper,  1971),  does  indeed  indicate  that  the  coupling  factor  for  peak  dis¬ 
placements  could  be  about  one- third  to  two- thirds  the  coupling  factor  for 
peak  particle  velocities.  Thus,  it  is  felt  that  it  would  be  conservative 
to  use  a  coupling  factor  of  about  3  percent  for  estimating  peak  radial  dis¬ 
placements  below  a  large  yield  surface  burst  on  rock  Thus,  for  5  percent 
coupling,  the  best  estimate  for  direct  induced  radiai  accelerations,  par¬ 
ticle  velocities,  stresses  and  strains  directly  below  a  surface  burst  on 
rock  are  given  by: 
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If  a  3  percent  coupling  factor  is  used  for  estimating  peak  radial 
displacements,  then  Eq.  35  may  be  altered  to  apply  to  the  surface  burst 
case  to  yield 
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Equations  36,  37,  38,  39  and  40  are.  applicable  for  estimating  radial  motions 
along  a  vertical  radius  extending  below  the  center  of  the  crater  of  a  surface 
burst,  i.e.,  8  =  0°  Fig.  5.  For  motion  predictions  along  radii  for  e  <  70°  , 
Fig.  5,  theoretical  studies  (Cooper,  Brode,  and  Leigh,  1971)  indicate  that 
Eqs.  36,  37  and  38  can  be  multiplied  by  a  correction  factor  F0  to  give 
reduced  motions  along  radii  at  6>0°.  The  value  of  F0  suggested  by  Cooper, 
Brode,  and  Leigh  (1971)  is  given  by 

FQ  «  (cosO  +  o.l  sin  0)  (41 ) 

where  0  is  defined  in  Fig.  5.  Since  Eq.  41  is  based  on  highly  complex  com¬ 
putational  models  (Trulio,  1970),  and  the  results -for  values  of  9  greater 
than  60°  are  very  sensitive  to  the  unloading  properties  assumed  for  the 
medium,  it  is  felt  at  this  time  that  a  value  of  F^  should  be  used  for  design 
which  is  more  conservative  than  the  value  given  by  Eq.  41.  The  following 
correction  value  is  suggested. 


0.5  (1  +  cos  o) 


(42) 
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li 


Equation  42  can  be  applied  to  Eq.  40  for  displacements  for  values  of  0 
up  to  65°.  but  cannot  be  used  for  values  of  0  >  65°  because  particle  dis¬ 
placements  at  shallow  depths  from  cratering  effects  are  not  predictable  by 
using  the  seismic  velocity  as  the  pertinent  property  of  the  medium.  The 
prediction  of  crater  induced  motions  at  shallow  depths  has  been  treated 
elsewhere  (Cooper,  Brode  and  Leigh,  1971).  It  should  also  be  pointed  out 
that  Eqs.  36,  37,  38,  39  and  40  should  not  be  used  for  media  with  seismic 
velocities  less  than  about  900U  ft/sec  because  the  tuff  data  (c  -  6000 
fps)  shown  in  Figs.  3  and  4  does  not  correlate  with  the  data  upon  which 
these  equations  are  based. 

It  should  also  be  noted  that,  at  this  time,  the  equations  developed 
herein  should  not  be  modified  for  use  in  design  cases  where  the  rock  medium 
in  question  is  saturated  and  below  the  water  table.  Much  of  the  data  in¬ 
cluded  in  this  analysis  has  been  from  rock  masses  which  were  saturated  and 
the  data  already  inherently  includes  the  effects  of  porepressures  if  they 
are  significant.  This  is  especially  true  of  the  Long  Shot  data  obtained  on 
Amchitka  Island.  Field  information  also  indicates  that  the  tuff  and  granite 
formations  at  Nevada  Test  Site  are  essentially  saturated  at  the  depths  at 
which  information  has  been  obtained  and  the  sandstone  formation  from  which 
the  Gas  Buggy  data  ware  ibtained  was  also  saturated. 
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TABLE  I 


Variables  Considered  in  Dimensional  Analysis 
of  Explosion  Phenomena 


Variable  Symbol 

(Independent  Variables) 

Energy  Released  by  Explosion  W 

Distance  from  Explosion  (Range)  R 

Seismic  Velocity  of  Rock  Mass  e 

Density  of  the  Rock  Mass  p 

Time  t 


(Dependent  Variables) 
Displacement  <$ 
Particle  Velocity  v 
Acceleration  a 


Dimension 
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L 
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TABLE  II 


Events  from  Which  Ground  Motion  Data  Were 
Obtained  and  Scaled 
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All  data  and  distances  scaled  by  cube  root  scaling  to  1  Kt. 
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